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Abstract 

We look for similarity transformations which yield mappings between different one- 
dimensional reaction-diffusion processes. In this way results obtained for special systems 
can be generalized to equivalent reaction-diffusion models. The coagulation (A + A — > A) or 
the annihilation (A + A — > 0) models can be mapped onto systems in which both processes 
are allowed. With the help of the coagulation-decoagulation model results for some death- 
decoagulation and annihilation-creation systems are given. We also find a reaction-diffusion 
system which is equivalent to the two species annihilation model (A + B —>■ 0). 

Besides we present numerical results of Monte Carlo simulations. An accurate description 
of the effects of the reaction rates on the concentration in one-species diffusion-annihilation 
model is made. The asymptotic behavior of the concentration in the two species annihilation 
system (A + B — > 0) with symmetric initial conditions is studied. 
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1 Introduction 



Diffusion controlled reactions in one dimensional systems have attracted much interest in 
the last years. They are non-equilibrium statistical systems which exhibit the property of 
self-organization. Mean field rate equations fail to reproduce the dynamics of these models 
which are characterized by non-trivial correlations. Theoretical descriptions of such systems 
have to take into account local fluctuations in the particle density. 

Exact results have been obtained for some one-dimensional models in the continuum from 
diffusion like equations ]I] and on a one-dimensional lattice 0. |3|. (|. In the latter case the 
time evolution of the system is determined by a master equation || . It is useful to rewrite it 
clS db Euclidian Schrodinger equation M. 

Interesting experimental measurements have been reported recently || [F| |[ [|. With the 
help of these experiments the exciton propagation is studied. The conclusion drawn from 
the measurements is that this propagation takes place in one dimension along structures like 
molecular "wires" and chains. The luminescence decay at early times is algebraic: 

c(t) oc i (1.1) 

with a = 1/2. This observation yields the conjecture that exciton-exciton annihilation 
processes take place |J. However, it is not clear whether exciton-exciton coagulation (fusion) 
reactions happen. In both the annihilation and coagulation chemical models the long-time 
behavior of the concentration obeys (|1 . 1| ) . Moreover, we will prove that this asymptotic 
behavior of the concentration is common to all models in which both reactions take place 
( |3.10| ). We call them annihilation-coagulation systems. The coefficient of the leading term 
depends on the on the ratio of the rates of these two processes. 

We show in this paper that all exciton reaction rates can be determined from experimental 
data. In our study we use a numerical and an analytical approach: 



First we want to determine the next to leading term in the long-time expansion of the 
concentration in the coagulation model. Three results are given in the literature concerning 
this regime. An approximative equation for the particle concentration derived from empty 
interval probabilities is obtained in [JITJ . A approximation scheme applied to an diffusion- type 
equations is used in |p |. These two approaches lead to almost the same conjecture: 



c(t) 




7 +0(t~ 



X2) 



where R — 1 in and R = ir/4 ~ 0.786 in [ 1Q[ . Here 1 — A' is the coagulation rate 



The time scale in which the diffusion rate is equal to one is chosen. The coefficient of the 
leading term was exactly computed for A' = [|l^]. In this case the next to leading term is 
of order t 3//2 [Q and the operator which appears in the master equation is the Hamiltonian 
of a supersymmetric system Q. 

We verify the conjecture ( |1.2|) by determining the value of the coefficient of the next to 
leading term in the A' ^ case. It can be computed by solving the Bethe Ansatz equations for 
the XX Z chain in an external field ||. However this computation is a mathematical challenge 
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because one needs to know both the spectra and the eigenvectors of the Hamiltonian. Here 
we investigate this problem with the help of Monte Carlo simulations. The numerical data 
presented in Section |] strongly support the R = 1 conjecture. 

Secondly, we want to find all chemical models which can be described with the help of 
( |1 . 2|) and other results concerning the coagulation model. We do this by looking for similarity 
transformations which map the coagulation model onto other react ion- diffusion systems. In 



12] , |l3j , |Tj] it was shown that the one-species annihilation and coagulation systems are con- 
nected. A general proof is given within the framework of the "Hamiltonian" formalism. The 
operators which appear in the master equations describing the two models are the Hamilto- 
nians of quantum chains || and are equivalent. We say that two reaction-diffusion systems 
are equivalent when the corresponding Hamiltonians can be obtained one from the other 
through a similarity transformation @]. In this paper we apply this method and look for fur- 
ther transformation which connect chemical models and lead to linear relations among their 
observables. Thus one can extend results concerning one reaction-diffusion systems to all its 
equivalent ones. 

We prove that the one-species coagulation-annihilation models can be divided in equiva- 
lence classes with respect to one parameter. Each class contains also a pure coagulation and 
a pure annihilation model. 

Another reaction-diffusion system which has attracted much attention is the two species 
annihilation model (A + B — > + 0). No exact results are available for this system. Scaling 
considerations indicate that, for equal initial site occupation probabilities of the two species, 
the concentration decay is algebraic with a time exponent equal to 1/4 In |16[ it was 



conjectured that the asymptotic behavior of the concentration depends on its initial value: 



c A (t) = c B (t) = K ^ / f ) + 0(r x ) if eu(0) = c B (0) = c(0), (1.3) 



where K = (27r)~ 3 / 4 ~ 0.252. This result is also obtained in |T7| from a renormalization 
group analysis for spatial dimensions higher then two. A general closure scheme for truncating 
the hierarchies of the joint density function equations leads to the same formula but with 
K = (327r)~ 1 / 4 ~ 0.316 PS|. A value of K — 0.28 has been numerically determined (16 



We report results of Monte Carlo simulations which allow us to obtain a better estimate for 
the constant K. It is in good agreement with the values conjectured in [TBI, IT? . 



We also look for the chemical models equivalent to the two species annihilation reaction- 
diffusion system. 

The paper is structured as follows: 

In section ^| a brief review of the "Hamiltonian" formalism is given. A general condition 
satisfied by the transformations which map reaction-diffusion system among themselves is 
derived. 

In section |3] we look for chemical models equivalent to systems previously studied. We 
concentrate on systems for which the particle concentration was computed. 

Two models which the procedure is applied to were mentioned before: the one species 
coagulation and the two species annihilation react ion- diffusion systems. Another model which 
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is exactly solved in some special case is the coagulation-decoagulation one f4|. We find 
two reaction-diffusion systems which are equivalent to it: the annihilation-creation and the 
death-decoagulation models. Thus we can give exact solutions for the particle number in a 
probabilistic cellular automaton in which decoagulation and death processes are permitted, 
but particles do not diffuse. 

In the last section new results obtained from Monte Carlo simulations are presented. 



2 The Master Equation and the Similarity Transforma- 
tion 

We begin with a brief review of the " Hamiltonian" description of reaction-diffusion sys- 
tems. We will use the same notations and conventions as in |4], |19[ where a comprehensive 
description of this formalism is given. 

First one has to construct the configuration space. Hard core interactions among particles 
moving on the one-dimensional lattice are considered. This also means that multiple occupa- 
tion of sites is not allowed. If the system is populated with one species of particles, each site 
can be either empty or occupied by one particle. If there are two species of particles each site 
can be in one of three states. The number of states is denoted by N and the lattice length 
by L. To each site % we attach a iV-dimensional vector space Vi. In the two state models 
(N = 2) the vector (1, 0) corresponds to a vacancy while the vector (0, 1) corresponds 
to a occupation of the site by an A particle. In the three state models (N = 3) the vectors 
( 1, 0, ), ( 0, 1, ) and (0, 0, 1 ) correspond to an empty site and to occupation of the site 
by a particle A or B respectively. 

The configuration space has the structure of an L-fold tensor product V — V\ ® • • • <8> Vl- 
Each of the N L vectors of the tensor product base corresponds to a possible configuration 
of the lattice. A natural convention used to denote these vectors is obtained by associating 
to each site a variable $ which takes integer values between and N — 1. In the two state 
models we use $ = to index a vacancy while = 1 corresponds to the presence of an A 
particle. In addition, (3i = 2 is used to describe a site that is occupied by a B particlei, in the 
three state models. In this 'spin basis' || 

\{(3}} = \(3 1} ...,(3 L ) (2.1) 

we define the ket vector 

\p) = E p ({Pht)\{P}) (2-2) 

W 

which describes the state of the system. P({f3};t) stands for the realization probability of 
the configuration {/3} = ...,/?£,} at time t. 

The state of a system can also be described by the set of empty intervals probabilities 
2H[ ] . They are defined as the probabilities to find one or more sequences of consecutive vacant 



sites. This description sometimes simplifies exact computations. 

The dynamics of reaction-diffusion systems are determined by the rates of the allowed 
processes. We consider only reactions which can take place when two hopping particles 
collide. So the set of parameters is given by the probabilities that a state (a,/5) on two 
adjacent sites will change into the state (7, 5) after one unit of time, denoted by 
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T*'g with the convention T*'p = 0. (2.3) 

Here a, [3, 7 and 5 take values between and N — 1. The sum of all reaction rates which 
modify the state of two neighboring sites is denoted by 

JV-l 

r^= £r#. (2-4) 

7,5=0 

All rates are non-negative, constant and real. We discuss only the left-right symmetric case. 
That is all the rates satisfy the condition 

r# = i#. (2.5) 

For one species systems a list of reactions and corresponding rates is given below: 

diffusion A + -> + A rate rj? 

annihilation A + A -> + rate 

creation + -> A + A rate r°° 

coagulation A + A A + rate 

decoagulation A + -> A + A rate T}? 

death A + -> + rate rjg 

birth + ^A + rate 

The master equation || describing the time evolution of the probability distribution 
P({/3}; £) can be written in the form of a Euclidian Schrodinger equation || 

||P) = -H\P) . (2.6) 

i7 operates on V . We consider periodic boundaries. Due to the fact that the reactions take 
place only between particles placed on two neighboring sites this Hamiltonian can be written 
as a sum: 

H = £ #i , (2.7) 
where Hi acts locally on V,® V^+i. This operator is defined by: 

iV-l 1 

#i = £ [r Q) ^r^- E r 5rt]. (2-8) 

o,/3=0 7,<5=0 

where E kl are N x N matrices with the entries (E H ) nm = 5k^i,m- This Hamiltonian is not 
necessarily hermitian. It ensures the conservation of probability and thus has a bra ground 
state 

(0\H = (2.9) 

which is the sum of all basis vectors 

(0| = £({/3}|- (2.10) 
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In other words : if we consider Hi as a matrix, equation ( |2.8| ) with ( |2.4| ) states that the sum 
of the entries from each column is 0. H has the same property (|2.9| ). 
The formal solution of ( |2.6|) is: 

|P(t)> = exp(-Pft)|P > (2-11) 
where the initial state is denoted by \Pq) ||. 

The expectation value of an observable X is given by the matrix element 

<X>(\P ),t)=j:X({(3})P({(3};t) = (0\X\P(t)}. (2.12) 

w 

The advantage of this Hamiltonian formalism is that we can introduce similarity transfor- 
mations. Our interest is to find those B for which the transformed Hamiltonian 

L 

H = Y, H i = BHB- 1 (2.13) 

i=l 

describes also the dynamic of a reaction-diffusion system. (The quantities in the transformed 
model are denoted with a *.) Therefore we have to concentrate on local transformations: 

B = b(g)b(g) ■ ■ ■ (g)b (2.14) 

defined by the L-fold tensor product of a real N x N matrix b. The mapping can easily 
be extended to the computation of expectation values. They can be computed M from the 
corresponding expectation values in the 'original' model: 

<X>(|P ),t) = (0\Xexp{-Ht)\P o ) = 



(0|X J Bexp(-Pt)p- 1 |P ) =< XB > (B _1 |Po>,*) (2.15) 

Some restrictions have to be imposed on b. First, the sum of the entries from each column 
of Hi must be 0: 

(0\BHi = 0. (2.16) 

A general solution of (|2.16| ) is: 

(0|P = p(0| (2.17) 

Here p is real. We have checked explicitly for every model considered in this paper that 
transformations obeying ( |2.16| ) but not ( |2.17|) do not lead to new mappings. Using (|2.10|) and 
( [2.14 ) one gets the restriction: 

N-l 

£&y = pi. (2.18) 

We can eliminate a overall scaling parameter by choosing p = 1. Thus the transformations 
are represented by real N x N matrices with the property that the sum of the entries on each 
column is equal to 1. 



5 



Not all these N x (N — 1) - parameter transformations lead to another reaction-diffusion 
system. The matrix elements of the two sites ( L = 2) Hamiltonian H have to satisfy the 
following sign conditions: 

Hmm >0 m = 0,l...iV 2 -l 

H mn <0 m,n = 0, 1...N 2 -!, m^n (2.19) 

The remaining transformations map reaction-diffusion system one onto the other. The pro- 
gram REDUCE was used to compute the transformed Hamiltonians. 



3 Mappings between reaction-diffusion systems 

We will start with the Hamiltonians of those chemical models for which the behavior of 
the particle concentration is known and look for local transformations which map them onto 
other reaction-diffusion systems. 

There are always the N\ — 1 permutation transformations which interchange (if N = 3 
different types of) particles and vacancies. We will treat everything "modulo" these permu- 
tations. 

As we mentioned before, we restrict ourselves to the left right symmetric case. The 
similarity transformation preserves this property. The generalization to the non-symmetric 
case is straightforward. 

In the three state model we take equal diffusion rates for the two kinds of particles (T°q = 
T^o). One parameter of the problem can be eliminated by choosing the time scale in which 
the diffusion rates are equal to 1. 



3.1 Two-state models 



Exact results are available for the one-species annihilation, coagulation and coagulation- 
decoagulation models. For the last two models the system of linear differential equations for 
the empty interval probabilities was written and solved. These equations decouple from those 
corresponding to many intervals probabilities [19|, 2T| . In the transformed models identical 



systems of equations can be written but for other observables, namely the probabilities of 
having interpolating sequences of vacancies and particles on consecutive sites . 
The general form of the local transformation is given by ( |2.18| ): 




(3.1) 

We will start with the annihilation model and not with the equivalent coagulation one. 
For simplicity we will use the notation a , instead of Tqq, for the annihilation rate and / for 
the coagulation rate (T\l). 
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The annihilation model. In this case particles diffuse and disappear pairwise when two 
of them try to occupy the same site. We denote the corresponding Hamiltonian with H ann 



and introduce, as in 0, ||, |25|, the parameter A' defined by: 



A' = l-~. (3.2) 

If a < 2, A' is equal to the probability that no reaction takes place at an two particles 
encounter. If the parameter A' is zero exact solutions can be obtained because H &nm can be 
expressed in terms of free fermions. In the time dependence of the particle concentration 
is derived in the thermodynamical limit with a full lattice as initial configuration. The finite 
chain is treated exactly in [[| and the two point correlation functions for the infinite chain are 
obtained. The particle concentration for random homogeneous initial conditions is deduced 
in [J- 

It can easily be seen that the transformed Hamiltonian H will only fulfill the conditions 
( |2.19D if at least one of the entries of b is nonzero. We take a = 0. The case (3 = can be 



obtained by applying a permutation transformation on H. 

It is a well known fact |12| , |13| , |l4jl that the annihilation and coagulation models are 
equivalent. We reobtain this result by choosing (3 = — 1 . In this case H is identical with the 
Hamiltonian of the coagulation model |4[|. 

This is, however, only one of the possible choices for the transformation parameter. For 
— 1 < (3 < we get the Hamiltonian of a system in which apart from diffusion with the same 
rate = 1, the following processes are permitted: 

• annihilation with rate 

f £ = a = 2^(1 - A') 



1-/5' 



coagulation with rate 



for = / = 2^-7(1 - A') 



The matrix b and A' can be written as functions of a and /: 

A' = 1 - / - a/2 (3.3) 
or equivalently a = a + 2f. One gets the following result: 

H(aJ)=BH ann B- 1 

where 

S|) (3-4) 

\ a+f I 

and with the parameter A' ( |3.3| ). 

We have obtained a one-parameter group of similarity transformations which act between 
annihilation-coagulation models ((a, /) — > (a, /)) corresponding to the same value of A' (i.e. 
a/2 + / = a/2 + / ). It is straightforward to write the transformation which connects the 
coagulation model (a = 0) and one described by H(a, /). 
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Some explicit results concerning the particle concentration can be given now by applying 
transformations Q2.15 ) on formulas obtained for the annihilation (/ = 0) ( p.4|) or coagulation 
(a = 0) models. 

Uncorrelated initial conditions are considered: at time t — every site is occupied with 
the same probability p. In this case: 



B- l \Po) 



'I - 



a+2f 



p\ 



_ a+f 
a+2f 

a +/ p 



p\ 



a+2f r i a+2f L L 

which corresponds to an initial site occupation probability in the equivalent annihilation 
fP- 



model equal to 



a+2f 



The operator of the occupation number of the site i is given by: 




1 



(3.5) 



The expectation value of the product of k such operators in the coagulation-annihilation 
model is (fZTJ): 

k 



< n il ---n ik >{p,t) 



a + V 



<n h --- n ik > E 



a + f 
,0 + 2/ 



P,t 



(3.6) 



We now drop the superscript. 

The exactly solvable case A' = l — / — f = corresponds to the physical situation that 
whenever two particle meet at least one of them disappears from the system. We mentioned in 
the introduction that the operators which appears in the master equation of these models are 
the Hamiltonians of supersymmetric systems. This reaction-diffusion models are mentioned 
in p2| in connection with the Glauber dynamic of the q state Potts model. The results given 
in M for / = 1 or a = 2 can be generalized to these annihilation-coagulation models. The 
particle concentration on a finite chain is: 

1 - [1 - (2 - /)# 



c(t,p, L) 



{2-f)L 

1 L - x 
(2=7)1 £ 



-l) fc [l-(2-/)# 1 - (-l) fe [l - (2 - f)p\ 



l+[l-(2-/)p]2 _ nk 
2[l-(2-/)p] tUb L 



1 - cos f 



irk 



x sin 

The behavior of the concentration 

L — > oo, t 



T ,exp 



irk 

At 1 - cos — 



L 

in the finite-size scaling limit ||: 
At 



oo with z 



L 2 



fixed 



(3.7) 



(3.8) 



can be determined Exact expressions can be written for the scaling and the first correction 
functions in terms of Jacobi theta functions. We get the following scaling relation: 



+5" 



2 ±a„ fn 



(2 - f)pj dz ° 3 I ' 2 J 



(3.9) 



Another interesting result concerns the large time behavior of the particle concentration 
in the thermodynamic limit. In the A' ^ case we start from the conjectured expression 
( |PD . Applying (pTToD we get: 



c(t) 



a + 2f 
a + f 



II R A' 1 
8rt + 2tt 1 - A' t 



+ o{r 



(3.10) 



The coefficient of the leading term depends only on the ratio of the annihilation and 
the coagulation rates. The next to leading term (~ - ) is dependent of A'. Both of them 
are independent on the initial concentration. We will give strong numerical evidence which 
supports the conjecture presented in [|ITJ (i.e. R = 1) in Section (J|). 

This asymptotic expression is an exact result in the supersymmetric (A' = 0) case 0. 
The next to leading term is of order t~ 3 ^ 2 . For uncorrelated homogeneous initial conditions, 
in the z — > limit M of fl3.9|) one gets: 



c(t,p) 




(2 - f)p 



1 



+ o(t- 5/2 ) 



(3.11) 



The coagulation-decoagulation model. In this system, apart from diffusion, coagulation 
and its backward reaction are permitted. We will use, as in the variable r\ defined through 
the decoagulation rate: rf = T\1 + 1. 

Exact solutions of this model are available in the continuum limit [0] and for the finite 
lattice in the case in which the coagulation and diffusion rates are equal (/ = 1) Q. In the 
latter case the steady state concentration is c s = 1 — r]~ 2 , independently of L and p. 

We found that the coagulation-decoagulation model is equivalent with two other reaction- 
diffusion systems. Only two transformations, which correspond to special choices for a and 
j3 in ( ft.l| ), connect this model with others. 

a) One of them is: 

6= (oT + 4) (3 ' 12) 

If / < 1) the transformed Hamiltonian H corresponds to a reaction-diffusion system with 
the following processes: 

• diffusion with rate = 1 — f 

• decoagulation with rate r\1 = f + rf — 1 

• death with rate rJo = / 
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If f = 1 there is no diffusion in the new model. The motion of the particles is realized 
by successive processes of decoagulation (with rate rf > 1 ) and death (with rate 1) on 
neighboring pairs of lattice sites. In this case we can use the results from [Q . For the finite 
chain and homogeneous random initial conditions the particle concentration is: 



2\lL 



c(t,p,L) 



i-[i- P (i-?r 2 )] 



v 



-2L 



L-l 



1 + (-l) k+ V[l - K 1 - V~ 2 )] L 



l+77 2 [l-p(l-rj- 2 
2ri[l-p(l-rr 2 ) 



COS ^ 



x sin 



^)exp(-A^)}. 



where A° is: 



A° k = rj\2{rj + rj~ l ) - 4 cos 
The full lattice is a stationary state: 

c(0,l,L) = c(t,l,L) = 



1. 



-l) fc+ V[l -p(l -r/- 2 )] Ll 



r? 2 +l 
2 V 



(3.13) 
(3.14) 



Interchanging A and a system with a death rate ?] 2 greater than the decoagulation rate 
1 and zero steady state concentration is obtained. For this choice of rates ( |3.13| ) gives the 
concentration of vacancies if p is replaced by 1 — p. 



b) The other transformation which gives a mapping of the coagulation-decoagulation 
model is: 



/ l 



V 5 



1 + 



^ + 1-1 



(3.15) 



One obtains an annihilation-creation reaction-diffusion system with rates: 



diffusion fJi 



v 2 +i 



annihilation f = a 



+ 1 + 1 



./" 2 



creation T^? 



n 2 



./ 2 



The annihilation rate is greater than the creation rate. The time unit can be adjusted in 
such a way that the new diffusion rate Tq° becomes equal to 1. 

In the case / = 1 the sum of the annihilation and creation rates is equal to 2. We can 
redefine the variable rj: 

2 

T] = 1 + 



2-o 



1 
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where y = 1/2 . The concentration for the finite chain is: 



c{t,p,L) 



1 - r] (77 - r] l ) [1 - (1 - 2p) L ?7 L ] 



2L 



L-l 

E 

fc=i 



2(1- V - 2L ) 
l + (-l) fc+1 (l-2p) L 



1 + 



l-2p) 



l+(l~2p) 2 
2(l-2p) 



COS 



nk 



V 2 +l 
2r? 



COS • 



x snr f ^- ) exp ( A"/) 



where A° is: 



2r] 

T] 2 + l 



2(rj + T] ) — 4 cos 



T 



(3.16) 



(3.17) 



The steady-state concentration is c s = (1 — rf v )j2 < 0.5. 

After a permutation of the two lines in (|3.15|) a system with interchanged annihilation and 
creation rates is obtained. The redefinition of rj in the exactly solved case is the same as the 
one given above but with y = —1/2. Formula (|3.16| ) gives the vacancies concentration if p is 
again replaced with 1 — p. The steady-state concentration is c s — (1 + r]~ 1 )/2 > 0.5. 

In the limit 7] —>■ 00 the reaction-diffusion system in which all three rates are equal is ob- 
tained. The operator which appears in the corresponding master equation is the Hamiltonian 
of the Ising model ||. 



3.2 Three state models 

The starting point in this section is the two-species annihilation model. The particles 
diffuse and react when an A and a B try to occupy the same site: 

A + B — > + 0. 

We denote the annihilation rate with a ( Tqq = a ), as in the preceding section. All diffusion 
rates are equal to 1 and = modeling a hard-core interactions between particles. 
The two sites Hamiltonian for this model is given by: 



( 














—a 





—a 


\ 





1 





-1 























1 











-1 











-1 





1 



























































a 

















-1 











1 





























a 





V 























/ 



(3.18) 



inthe{|00), |01), |02), • ■ ■ |22)} base. 

A general transformation has the form ( |2.18 ): 
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/ 1-a-p 



V 



a 
P 




(3.19) 



We have to insure that the transformed Hamiltonian satisfies the condition ( 2.19Q . The first, 
fifth and ninth columns of Hi are identical up to a factor. So are the line one, five and nine. 
The matrix elements of H can be interpreted in terms of reaction rates only if the entries 
from at least one of these lines or from its corresponding column are zero. The other two 
lines i,j G {1,5,9} have at most two non-zero entries in the two corresponding columns 
By imposing this conditions we are left with one transformation: 




(3.20) 



which leads to a new react ion- diffusion system if a ,a and j3 are greater than one and if 
or x + P^ 1 > 1 . The transformed two sites Hamiltonian is: 



H = 2 x 



1 














-a(-l + a- 1 + [3- 1 ) 





-a(-l + a -1 + P' 1 


o \ 





1 





-1 





-(1-a- 1 ) 





-(a - 1)(1 - a" 1 ) 











1 








_( G _l)(l_ r i) 


-1 










-1 





1 





-{a- 1)(1 - a" 1 ) 





-a-*- 1 ) 















































a 

















-1 








-(i-r 1 ) 


1 


_(a _!)(!_ 0-1) 


























a 





V o 























o / 



(3.21) 

H has the same diffusion part and is left-right symmetric. All processes which reduce 
the number of particles start from a AB pair. Apart from annihilation we have 'two species 



coagulation' reactions: 












• A + B - 


^0 + 


with rate 


fl2 
1 00 


= a (a' 1 + fi- 


-l _ 


1) 


• A + B - 


^0 + A 


with rate 


f 12 
1 01 


= I -a' 1 






• A + B ~ 


-> A + 


with rate 


f 12 
1 10 


= (a-l)(l- 


- a~ 


') 


• A + B - 


-> 5 + 


with rate 


pl2 
1 20 


= l-P' 1 






• A + B - 


^0 + 5 


with rate 


fl2 
1 02 


= (o-l)(l- 


-P~ 


- 1 ) 



Note that the diagonal terms of the Hamiltonian are invariant under the transformation 
( Ti2 = T 12 = a). The coagulation rates depend on the type of the particle which survives 
the process and they are (a — 1) times smaller if the particle which survives is the one 
which jumps. There are only two coagulation rates, independent of the initial position of the 
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surviving particle if a = 2 . This case corresponds to the physical situation in which at any 
A B encounter at least one particle leaves the system. 

From fl2.15| ) and ( |1.3|) it is easy to see that the concentration will have an algebraic fall-off 



if: 



ca(0) _ a 
cb(0) ~ P' 
In this case the particle concentrations are: 



c A (t)=a* V tl/4 W (3.22) 



c B (t) = /3* \ K> = ^c A (t) (3.23) 

4 Monte Carlo Simulations 

Only a few chemical models have been solved exactly. For some others, quantitative 
estimates have been made with the help of approximation schemes. This is the reason why 
numerical methods have been extensively used in the study of reaction-diffusion systems. The 
dynamics of the particle density and spatial distribution can be determined with the help of 
Monte Carlo simulations |23], [24[. 

The Hamiltonian formalism opens new possibilities in this field. For small chains numerical 
data can be obtained with a high accuracy by using diagonalization techniques or simulations. 
The study of the finite-size scaling behavior (|3.8|) in the limit z — > permits the determination 



of the particle concentration for infinite chains at very large times [25]. 

In this section we present our numerical results concerning the asymptotic behavior of the 
particle concentration in the one and two species annihilation models. We simulate directly 
the thermodynamical limit by using large lattices. 

We consider chemical systems in which no reaction is permitted which creates particles 
on a pair of empty sites. This is why the so called 'direct' Monte Carlo Method ||23| , |26| can 



be used. The way we implement this method is described in detail in j25 . 

We will first present our results for the two state models and continue with the three state 
models. 



4.1 The two state annihilation-coagulation model 

The influence of the reaction rate in annihilation and coagulation models is a subject 
which has received a lot of attention in recent years (see [R], |TT| and references therein). We 
mentioned in the introduction the conjectured expression (|1.2|) for the concentration decay 
in the A' ^ coagulation model. This result is in agreement with the conclusions of a 
numerical study we carried out in our previous work. Using a finite-size scaling analysis we 
obtained the same value for the time exponents in the leading (—1/2) and the next to leading 



term (—1) and we found that the coefficient of the leading term is independent of A' [25 . 

As mentioned there, a better test of (|1.2p (or equivalently of ( |3.10|) ) can be made through 
a direct simulation of the thermodynamic limit. Results of such computations are presented 
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in pn|| . They are, in some cases, in semi-quantitative agreement with the conjecture presented 
there. We concentrate on the long time behavior and use smaller lattices as the one mentioned 
in [11]. Thus we are able to obtain very good statistic in a reasonable CPU time. 

We tested the conjecture ( |1.2|) in two different ways. Simulations of the annihilation model 
were performed on a lattice of length L = 2000 for times up to a value of t max = 1.5 x 10 4 . 

In the first set of simulations two values for the initial concentrations, 1.0 and 0.1, were 
used. A decrease of 3 respectively 2 orders of magnitude of the concentration enables us to 
reach the asymptotic regime. 300,000 runs were performed in order to insure a relative error 
of the particle concentration of less then 10~ 3 . In Figure 1 the following quantity (see (|3.10|) 
with / = 0): 

BW= 2.(c( t )-^)^ (4.1) 

is represented for A' = 0.875 and 0.75. We see that our data converges to the value R(t) = 1 
for which the straight line is drawn. 

In Figure 2 the same kind of data corresponding to A' = —7 is shown. The convergence 
is worse because the absolute value of the correction term fl3.10| ) is 8 times smaller than for 



A' = 0.875 so R(t) approaches faster the order of magnitude of the numerical errors. 

A second set of simulations of the annihilation model were performed in which we averaged 
only over 20,000 runs. We applied a \ 2 test in two steps to the data. A linear combination 
of time powers was presumed: 

y = K x t~ x + K x+l/2 t- x ~ 1 / 2 + K^t-*- 1 (4.2) 

The first fit was made for the concentration y = c(t) taking x = 1/2. The coefficient of 
the leading term was determined with a confidence level greater then 99%. The first three 
significant digits of K1/2 are identical with those of (87r) -1 / 2 ~ 0.1995. 

The next step is to determine the coefficient of the first correction. The leading term is 
subtracted from the numerical data: 

y = c(t) - — L 



and a value x — 1 is used in the fit ( |4.2| ). The confidence level was good (greater then 50%). 
The values for R are listed in table |. 

For most values of A' simulations of the same model with different initial concentrations 
were performed, as can be seen in table |. The best fits were obtained for the choice of 
reaction rates corresponding to a value of A' near 1. This is easy to understand because for 
such A's the contribution of the next to leading term in (|3.10|) is maximal. The simulations 



corresponding to A' = —7, 3/4, and 7/8 are the ones used to make the Figures 1 and 2. 
For this data the statistic was made by averaging over 300,000 runs. 

The conclusion we draw from the results of the two tests is that the conjecture presented 
in |llj, predicting a value of R — 1 in (|1.2|) is correct. 



This results are valid if the data we considered before are not affected by the finite size 
of the lattice on which the simulations were performed. We check this by performing the 
simulations described in the following. 
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A' 


Initial concen- 
tration r(0) 


Values of i? 

V LAi± IAVjU V.-' J- J- I, 

from 1VT data 


7 
— 1 


i nn 
l.UU 


i non -U n non 
l.uyu ± U.UZU 


7 


n in 
U.1U 


i nfin i n nm 
l.UDU ± U.U1U 


Q /O 
— 6/ A 


1 nn 
l.UU 


n oa n i n nin 
U.»4U ± U.U1U 


Q7 /OH 
— Z / /ZO 


i nn 
l.UU 


n oe;n -u n non 
U.yoU ± U.UZU 


1 

— 1 


i nn 
l.UU 


i nnn x n n<in 
l.UUU ± U.U4U 


1 

— 1 


n 9^ 


i i on -U n nsn 
l.loU ± U.UoU 


1 

— 1 


n i n 
U.lU 


i nRn -U n nsn 
l.UUU ± U.UoU 


— 3/7 


1 nn 

l.Wu 


1 m n + n n3n 

1.U1U 1 u.uuw 


-1/9 


l.nn 


n.94n ± n.n7n 


1/2 


l.nn 


1.060 ±0.040 


1/2 


n.25 


0.990 ± 0.040 


3/4 


l.nn 


1.004 ±0.005 


3/4 


n.in 


1.030 ±0.020 


7/8 


l.nn 


0.998 ±0.005 


7/8 


n.in 


1.005 ± 0.005 



Table I: Values of R (11.21) obtained from Monte Carlo simulations of the annihilation model. 



The two different conjectured values are 1 [11] and 0.786 10 



In the asymptotic regime the system contains a low number of particles which are spatially 
separated. A particle which was not involved in a reaction represents a random walker. This 
means that at large times, of the order of the square of lattice length, finite-size effects should 
be dominant. 

We obtain a more precise quantitative picture from simulations of the annihilation model 
for times up to t max = 10 6 . We use a lattice of the same length ( L — 2, 000 ) and averaged 
over 20,000 runs. Three values for the parameter A' were chosen: 7/8,0 and— 7. Each 
simulation started with a full lattice as initial configuration. 

For an infinite system the long time behavior of the concentration is algebraic . 

For a finite system the concentration decay is exponential. We are interested to study the 
transition between the two types of decay, i.e. the onset of finite size effects at large times. 

Figure 3 gives a log-log plot of the concentration. The straight line corresponds to c(t) = 
(87rt) -1//2 . We see that for times t > 10 5 the dashed lines strongly deviate from the straight 
line. 

In Figure 4 which is a plot of the logarithm of the concentration we can see that starting 
from the same value of time (t = 10 5 ) the curves fit nicely with another straight line described 
by: 

c(t) = l e xp(-2^t) . (4.3) 

The formula is obtained in the t — >• oo limit of ( p.7|) , for the A' = annihilation model 
(c=0). 

From Figures 3 and 4 we conclude that the onset of the finite-size effects takes place in a 
narrow vicinity of t = 10 5 . This value is considerably smaller than L 2 but much larger than 
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the t max we used. It is now clear that the results we presented in the first part of this section 
refer to the thermodynamic limit of the annihilation model. 

The slope of the line from Figure 4 is equal to the energy of the first excited state of 
-f^ann.- We can conclude that this energy has the value 4(1 — cos(7r/L)) independent of the 
parameter A'. This confirms previous results obtained from numerical diagonalization of this 



Hamiltonian [27 



4.2 Numerical results for the three state models 



The simulation of the thermodynamic limit of the two species annihilation reaction requires 
a very large amount of CPU time. 

This has two reasons. A reliable study of the long time decay of the concentration presumes 
simulations in which the total number of particles reduces by at least two orders of magnitude. 
The algebraic fall-off is much slower than in the A + A — > + case. This means that 
we have to perform the simulations up to much higher values of t max than we did in the one 
species case. 

On the other hand, if we start at t — with a random distribution, there will be local 
fluctuations of the particle densities. The decay of the initial fluctuations plays an essential 



role in the dynamics of the system fll6| . Their length scale extends in time. So one is forced 
to use longer lattices than in the one-species case where such domains do not appear. 

A lattice of length L = 10 5 was used to simulate three state models. In the case of the 
two species annihilation we were able to average only over 100 runs. Our aim is to see if the 
particle concentration obeys the algebraic decay law (|1.3|) for the A — B symmetric case and to 



compute K. Results from similar simulations were reported in ||T6[ which lead to K = 0.280. 
This is in equally good agreement with the two conjectured values 0.252 and 0.316. 
Figure 5 shows the curve of 

K{t) = c A (t) 4= (4.4) 
/c(0) 



for various choices of the annihilation rate and initial site occupation probabilities. The 
simulations were stopped at t max = 10 6 , a value hundred times larger as the one used in (16 



At very large times the curves converge showing that the asymptotic regime was reached. 
The average value is K = 0.247 ± 0.004 . 

Our result is important because it confirms the conjectured dependence of the concentra- 
tion on its initial value |16|]. It also suggest that the asymptotic behavior is universal in the 



sense that it is independent of the annihilation rate a. 

We cannot make any prediction concerning the exponent x in ( |1.3| ). This would require 
a better statistic in the determination of the asymptotic behavior of the concentration. This 
is an aim which we cannot achieve with the computer facilities which are currently at our 
disposal. 
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Conclusions 



We found a general condition satisfied by local transformations which realize mappings 
between reaction-diffusion systems. They are iV x (N — l)-parameter real matrices which 
have the sum of entries in each column equal to one. 

Some new results concerning two state models are derived: 

All annihilation-coagulation systems corresponding to the same parameter A' ( |3.3| ) are 
equivalent. The continuous similarity transformation which acts between these models is 
given. Results previously derived for the coagulation and annihilation models are generalized 
to the chemical systems in which both reactions are permitted. 

At this point we return to the exciton decay measurements mentioned in the introduction. 
Further experimental determinations are very desirable for two reasons: 

1) It would be interesting to have accurate measurements of the exciton hopping time 
which is related to their diffusion constant. The coefficient of the leading term ( ~ t -1 / 2 ) in 
the exciton concentration decay can then be extracted from the experimental data ||. Thus 
one can determinate the ratio of the exciton-exciton annihilation and coagulation rates ( p7T0|) . 



2)Lower errors in the experimental measurements of the luminescence decay would permit 
the computation of the next to leading term (its time exponent and coefficient) which appears 
in an expansion of the exciton concentration. 

If this term is of order t~ 3 / 2 (pj.ll[) the dynamic of exciton propagation and reactions 



is described by a system which is supersymmetric ||. This would be the first such system 
experimentally measured. The other possibility is that the next to leading term is of order 
t^ 1 ( |3.10| ). In both cases the determination of the leading and next to leading term and of 
the hopping time would permit the computation of the rates of the processes which take place 
between excitons. 

The two states coagulation-decoagulation, annihilation-creation and death-decoagulation 
models are equivalent. This enables us to give an exact expression for the particle number of 
a probabilistic cellular automaton in which modifications of the configuration occur only by 
starting from particle- vacancy pairs (|3.13|). The same observable is computed for the model 



in which the sum of the annihilation and creation probabilities is one (|3.16|) . 

A reaction-diffusion system equivalent to the two species annihilation ( A + B — > + ) 
is found. 

We have accurately determined the effect of the reaction rate on the long time decay of the 
concentration in the one-species annihilation model. Our results are in excellent agreement 



with a theoretical conjecture previously made in \TL \. 

We also present quantitative analysis of the asymptotic behavior of concentration in the 
three states annihilation model. They are also fully compatible with previous results. Al- 
though they are, as far as we know, the best available in the literature, they do not allow us to 
make any prediction concerning the effects of the reaction rate as we did for the A+A — > 0+0 
model. 
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List of Figures 



The Monte Carlo data are represented by interrupted lines. Error bars are given only for 
a few points. 

Fig. 1: Effects of the annihilation rate on the asymptotic behavior of the particle concentra- 
tion for the A + A — > + reaction - R(t) as defined by (|PD , for A' = 3/4 and 
7/8. The straight line is R(t) = 1. 

Fig. 2: Effects of the annihilation rate (2). R(t) as defined by (fOP , for A' = -7. 

Fig. 3: Onset of finite size effects - the algebraic decay. Double logarithmical plot of the 
concentration for the one-species annihilation reaction. 

Fig. 4: Onset of finite size effects - the exponential decay. Plot of the logarithm of the 
concentration for the one-species annihilation reaction. 

Fig. 5: Asymptotical behaviour of the particle concentration in the A + B — ► + reaction. 
Plot of K{t) as defined by (|4.4j) . The straight line corresponds to K(t) = (27r)~ 3 / 4 
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